Hole spectral functions of LaMnOs 
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By use of the orbital t-J model, we calculate the photoemission spectra of LaMnOa using the ex- 
act diagonalization technique, and interpret our numerics quite well in the orbital-polaron scenario 
where the scattering between holes and orbital excitations is treated within the self-consistent Born 
approximation. The quasiparticle bandwidth is found to be of the order of J and t in the purely 
Coulombic and Jahn-Teller phononic model, respectively. We suggest that angle-resolved photoe- 
mission spectroscopy experiments allow one to distinguish between the orbital-polaron scenario and 
the Jahn-Teller polaron scenario. 
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Orbital Physics is a key concept for recent intensive 
studies on transition- metal oxides j^] , especially on man- 
ganese oxides with perovskite structure i?i_x^xMn03 
{R = La, Pr, Nd, Sm and A = Ca, Sr, Ba) due to 
the discovery of colossal magnetoresistance in this class 
of materials. Mn^+ in /JMnOa has four d electrons of 
which three are put into the t2g orbitals and form an 
S — 3/2 localized spin, and the mobile one occupies one 
of the Cg orbitals {d^^-y^ and d322_^2). This Cg orbital 
degree of freedom may cause new phenomena through 
strong coupling with charge, spin, and lattice dynamics 
For example, alternating orbital order was observed 
in the ferromagnetic (FM) planes of A-type antiferro- 
magnetic (AF) LaMnOs ||^ and uniform x"^ — orbital 
order in A-type AF Pri/2Sri/2Mn03 [Q, while orbital hq- 
uid was proposed to exist in metallic FM Lai-^^Sr^jMnOa 
{x > 0.2) Q. These discoveries attract attention to the 
orbital correlation, dynamics and order-disorder transi- 
tion. An essential first step is to understand the motion 
of a hole in an orbital ordered system such as the FM 
planes of LaMnOa Q. We address it in this Letter. 

Orbital ordering in LaMnOa can be induced by ei- 
ther the Jahn-Teller (JT) lattice distortion originat- 
ing from the degeneracy of the eg orbitals or 
the intra-atomic Coulomb interaction in the Cg orbitals 
P, P^|JTllJl^ , [l3|] . In the latter mechanism, the Coulomb in- 
teraction eliminates doubly occupied sites and results in 
the so-called orbital t-J model, a purely electronic Hamil- 
tonian descrbing the orbital sector of physics of LaMnOa. 
The model was initially derived by Kugel and Khomskii 
[0 and has been recently studied within mean field the- 
ory by Ishihara et al. 1 1 1 1 and via a finite-temperature 
diagonalization method by Horschet a/. |l2|. The exper- 
imentally observed orbital ordering in LaMnOa |^ was 
reproduced by use of this model ]l^ . However, it can 
be also attributed to Jahn-Teller polarons In this 

Letter, we take the orbital t-J model as a starting point 
and discuss the JT effect later. 

We first calculate photoemission spectra for the or- 
bital t-J model using the exact diagonalization technique 



(ED). The experimental counterpart of this problem 
is angle-resolved photoemission spectroscopy (ARPES) 
measurements. ARPES experiments have not yet been 
reported in this doping regime. Thus, we perform here 
computational experiments to explore unbiased informa- 
tion. Then, we satisfactorily interpret the outcome of 
such experiments in an orbital-polaron picture where the 
scattering between holes and orbital excitations is treated 
within the self-consistent Born approximation. Further- 
more, we point out that although the JT lattice distortion 
leads to a large gap in the orbital excitation spectrum and 
thus resists the formation of the orbital-polaron, the hole 
can move almost freely through the orbital-flip process. 
The quasiparticle bandwidth is of the order of J (t) in the 
purely Coulombic (Jahn-Teller phononic) model, respec- 
tively. Our results indicate that ARPES may provide a 
possible approach to distinguish between the Coulombic 
scenario and the Jahn-Teller scenario. 
The orbital t-J model is [|To|.RIlJill 



-H- [T?T^ + 3T^T? T V3{T^T^ -f T-^Tf)], (1) 
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where d 



erator for the eg electron at orbital a. T( 



) is the constrained fermion op- 
id\^d.^^ - 

dl^d.^)/2 and Tj" = {dl^d + dl^d ■-^)/2 are orbital- 
pseudospin operators with | f) = dx2_y2 and | i) = 
d^z^-r^. The anisotropic transfer matrix elements are 



tV 



y 3^/4, til 



t/A, and il/ = ^i/ = TV3i/4, here the 



=F sign distinguishes hopping along the x and y directions. 
The orbital superexchange interaction J= t^/{U—J') 
with U ( J') being the interorbital Coulomb (exchange) 
integral [jlll,^. In LaMnOa, the realistic parameters are 
estimated from photoemission experiments jlljlj] to be 
t - 0.72 cV, [/ ~ 5 eV, J' - 2 eV, thus % G.24t. 

The photoemission spectrum (PES) (^kadka)'^ El] 
is calculated by using the standard Lanczqs algorithm 
|]l6| . There are two bands and we find ('^kj^c'ki)'^ = 



1 



('^k+QT'^k+QT)'^ ^^^^ Q Fig. 1 shows the PES 

{dl^^dj^^)^ of the ground state of the 16-site cluster at 
half filling. At any k, there is a well-defined quasi- 
particle pole (i.e., zero orbital- wave) at the low energy 
side which is well separated from a broad, incoherent, 
multiple-orbital-wave background extending to the full 
free-electron bandwidth. The bottom and the top of the 
quasiparticle (QP) band locate at (0,0) and (tt, tt), re- 
spectively. The QP bandwidth of the order of J is much 
more narrow than the free-electron bandwidth. Note that 
all momenta belonging to the noninteracting 2D Fermi 
surface cos -\- cos ky = are no longer degenerate de- 
spite their close proximity in energy. Correspondingly, 
the van Hove singularity in the density of states is weak- 
ened in the orbital t-J model. 

Since the mixed terms cx T?=T.' + T-^T? in the Hamil- 
tonian (|l|) could contribute only in higher order orbital- 
wave theory than the linear one, it is interesting to per- 
form similar calculations using the ED technique for the 
simplified Hamiltonian without these terms, which has 
been referred to as "the truncated Hamiltonian" by Brink 
et al. in their study of the purely orbital J model 
This could let us investigate the contribution of higher 
order orbital excitations to the QP behavior. Fig. 1 also 
illustrates the PES for the truncated Hamiltonian. The 
overall line shape of the PES remains almost unchanged 
but a slight shift to a higher energy. The prominent in- 
fluence of higher order orbital excitations is that they 
suppress the heights of the quasiparticle peaks and thus 
weaken the QP spectral weights. We therefore conclude 
that as far as the QP dispersion relation is concerned, 
higher order orbital excitations can be neglected. The 
resultant QP dispersions are displayed in Fig. 2. 

To get the physical insight of the above numerics, we 
perform an analytical calculation of the QP spectral func- 
tions in an orbital-polaron scenario. The J term induces 
the orbital Neel state, i.e., the alternation of orthogonal 
orbitals in the ground state. Here we analogize our prob- 
lem to that of a hole moving in the AF spin background, 
an essential problem in the field of high-temperature su- 
perconductivity jl^. The latter can be accurately un- 
derstood in the spin-polaron scenario where holes are de- 
scribed as spinless fermions (holons) and spins as hard- 
core bosons. The scattering between holes and spin-wave 
excitations is treated within the self-consistent Born ap- 
proximation (SCBA) 0. In the SCBA, higher order 
spin-wave excitations are also neglected, consistent with 
our ED results. Hence, we expect that our numerics may 
be interpreted in a similar way. 

In LaMnOs, alternatingly occupied orbitals on the two 
sublattices are oriented in the FM planes: (| t) + l 
and (I t) - I i))/%/2 This state is different from 

the alternating directional orbitals — and iy^ — r^, 
which have been naively expected (ll]. To obtain the 
Neel configuration iTi^Ti^+iTi^+aTi+s • • ■> = I itit ' ' •> m 



a new orbital basis, we perform a uniform rotation of or- 
bitals by 90° about the T^' axis: di| (di| - dix)/\/2, 
dii ^ (diT + du)/\/2,and thus T? -Tf, T? . 
Note that in the new orbital basis, the interorbital trans- 
fer matrix element is changed to (<[/ — ^ij^)/2. It is the 
imbalance between the dx2_y2 — (i^2_j,2 hopping and the 
d^^2_j,2 — d^^2_j,2 hopping that gives rise to the interor- 
bital hopping in the rotated basis and causes the bare 
hole dispersion as shown below. Orbital excitations can 
be described within the linear spin- wave theory by defin- 
ing T'l = al, Tf = OjOi — ^ on the J, sublattice, and 
Tf = by T? = ^-bjb. on the t sublattice Here a-, 

and 6j are hard-core boson operators. Considering the 
orbital Neel state as the vacuum state, we define holon 
operators hi (/;) similar to those in the spin-polaron sce- 
nario |l7j so that di| — hla-^, di| = /ij on the | sublat- 
tice and dji = f Vby dji = / j on the f sublattice. 

Introducing new fermion operators {/k, h^} in the mo- 
mentum space /k = (/k + /ik)/%/2, /ik = (/k - ft.k)/V2, 
we arrive at an effective orbital-polaron Hamiltonian 

i^cff = J2 l/k -hih^) + EK<«q + </3^/^q) 

k q 

+ E(/k/k-q- /ik/ik-q)(.9fq/3q+5£q«q) +H.C. (2) 
kq 

+ E(^k/k-q - flhi^-ML^^ + -^kq^q) + H.C., 
kq 

where the bare hole dispersion is — ~^7k- The holon- 
orbital-wave coupling functions are ~ --2=('-yj^w+ -|- 

7k-qO. Pkq ^ -^(^k< - ^k-qO' P^q = .9k,q+Q 

and g£q = Pk,q+Q with 7k = (cos/c^; + cosky)/2, -q^ = 
(cos kx — cos ky)/2, and 

«± = {[(Aq±i?q)/c.± + l]/2}i/2, (3) 
< = -sgn(±i?q){[(Aq ± Bq)/c.± - l]/2}i/2. (4) 

Here Aq = 3J and Bq =Jj^^/2. The ofq's (/3q's) are 
orbital-wave operators, flq = -^(Uq/3q+i'q /^Iq+Uq Oq-I- 

Vqalq), bq = ^(uJ/3q + V + /3lq~Uqaq~'yq a^^q), with 

dispersion lj^ — y/Aq{Aq ± 2i?q). 

The holon Green's functions G/_?i(k, w) = =F Ek — 
i;/,/i(k, w)]"^ are calculated within the SCBA @. Note 
that S/(k-HQ,tj) = S/,.(k,w) and Gf{k + Q,uj) = 
Gh(k, oj) with Q = (tt, tt). The self-energy is thus of form 

EHk, CO) = E '[(5kq)'G'/(k - q, - O 
q 

+ (P^q)'G/(k + Q-q,c. -..+ )], (5) 

where means to sum over the first Brillouin zone. 

The QP spectral functions ^/ (k, w) — — ^ Im G/ (k, w) 
are shown in Fig. 3 The QP dispersion = 



2 



ek + ^fO^,El) (see Fig. 



Z(k) = 



1 - 



2) and the spectral weights 
(see Table I). All of the 

SCBA results are in good agreement with the ED re- 
sults, especially with those for the truncated Hamiltonian 
as expected. Therefore, the orbital-polaron scenario may 
provide a valuable scheme for further works on the orbital 
dynamics. 

In the rest of this Letter, we discuss the Jahn- Teller 
effect on the quasiparticle band. Recently, Bala and Oles 
suggested that a large gap in the orbital excitation spec- 
trum induced by the JT lattice distortion might lead to 
a strong confinement of holes in lightly doped LaMnOa 
insulators |]l3| . This implies that the hole quasiparticle 
bandwidth is more narrow in the presence of JT phonons. 
We include the JT effect in the same way as Bala and 
Oles [0,0 but we get a different result. The Jahn- Teller 
interaction considered here is 1131 



(6) 



Here T? refers to the rotated basis. Ejt{(P,5x,Sz,u) — 
A[((5^-(5^)sin2<?!)-2V^ u cos 20] acts as a fictitious "mag- 
netic field" in which (jj is the tilting angle of pseudospins, 
and 6x, Sz, u characterize lattice distortions (in units 
of the lattice constant): 5x{5z) — uniform deforma- 
tion along the x and y (z) directions, and u — oxy- 
gen ionic displacement along Mn-O-Mn bond in the xy 
plane |19|. The distorted lattice energy per Mn ion is 
Ei{5x.5z.u) = K,{\51 + 2u2 + \5l) + K^iSl + \6l) 
where Ki{K2) is the nearest-neighbor Mn-0 (Mn-Mn) 
spring constant. To estimate the JT effect on the quasi- 
particle behavior, we may consider the case of = 
(i.e., no static distortions due to a tetragonal field) with- 
out loss of generality. In the classical ground state at 
= 0, EjT = —3X^/Ki. For the realistic parameters of 
LaMnOa: the spring constant Ki — 200 eV and the JT 
interaction parameter A ~ 6 eV [|7|,p^, -Ejt = —0.54 eV. 

In order to calculate in the presence of JT lattice 
distortion, all we have to modify the above derivation is 
to make the following replacement in (|^) and (^: Aq 
3 J - 2EjT = 3J + eX^/Ki. The JT interaction adds 
an Ising-like component to the excitations and induces a 
large gap in the orbital excitation spectrum. Thus, the 
JT effect stabilizes the orbital ordering. 

The SCBA results on the 20 x 20 lattice are summa- 
rized in Table 1. The large gap in the orbital excitation 
spectrum induced by the JT lattice distortion does not 
weaken the QP bandwidth. Instead, at A = 6 eV the 
QP bandwidth reaches 75 percent of the width {2t) of 
the bare hole dispersion £k- Although the large gap in 
the orbital excitation spectrum resists the formation of 
the orbital polaron, the hole can move almost freely via 
the orbital-flip process in the rotated orbital basis. For 
Aq ^ t because of large J or large A, Eq. (||) can be 



solved analytically in perturbation theory. Most of the 
spectral weight, 1 — 0{t'^/Aq), appears in the quasiparti- 
cle part of G/(k, oj) which behaves indeed a bare disper- 



sion E, 



f 



Fig. 4 shows the QP bandwidth W as a function of J. 
Without the JT interaction (A = eV), W ~ 2.2J scales 
with J in the range of 0.01 < J < 0.4, and approaches to 
the width of the bare hole dispersion (2t) for large J. On 
the other hand, with strong JT interaction (A = 6 eV), 
W is of the order of t in the whole range of J. These 
results reminisce those for one hole moving in the spin 
t-t'-J model (the t-t'-Jz model with large J^), respec- 
tively, where the transfer to next nearest neighbors (t') 
provides the bare hole dispersion ||l^,|l6| . Therefore, the 
QP bandwidth for realistic J is of the order of t (J) in 
the presence (absence) of strong Jahn- Teller interaction. 
Our results can be tested by future ARPES experiments. 
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TABLE 1. Bandwidth W and spectral weights Z{k) as a 
function of the Jahn- Teller interaction parameter \. J — 0.3i. 
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FIG. 1. The PES (dl^dt^^)^ for the orbital t-J model on 
the 4x4 cluster using the ED technique for the full Hamil- 
tonian (solid lines) and the truncated Hamiltonian (dashed 
lines). Here J = 0.3t. 
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FIG. 2. Hole quasiparticle dispersion for the orbital t-J 
model with J = 0.3t: The ED results for the full Hamiltonian 

(solid squares) and the truncated Hamiltonian (solid circles); 
The SCBA results on the 4x4 cluster (open circles) and 
the 20 X 20 one (open squares) . Solid lines are the fits using 
eI= ao + oi7k + 02 cos kx cos ky + 0372k- 
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FIG. 4. The quasiparticle bandwidth W as a, function of J 
at different Jahn- Teller interaction parameters: A = (solid 
squares) and A = 6 eV (open circles). Calculations are per- 
formed on the 16 x 16 cluster within the SCBA. 
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FIG. 3. The PES A/(k, w) for the orbital t-J model on the 
4x4 cluster calculated within the SCBA. Here J = 0.3t. 
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